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ABSTRACT 

We perform high-resolution hydrodynamic simulations of a Milky Way-mass galaxy in a fully cosmo- 
logical setting using the adaptive mesh refinement code, Enzo, and study the kinematics of gas in the 
simulated galactic halo. We find that the gas inflow occurs mostly along filamentary structures in the 
halo. The warm-hot (10^ K < T < 10^ K) and hot {T > 10^ K) ionized gases are found to dominate the 
overall mass accretion in the system (with M = 3-5 Mq yr""'^) over a large range of distances, extending 
from the virial radius to the vicinity of the disk. Most of the inflowing gas (by mass) does not cool, and 
the small fraction that manages to cool does so primarily close to the galaxy {R < 20 kpc), perhaps 
comprising the neutral gas that may be detectable as, e.g., high- velocity clouds. The neutral clouds 
are embedded within larger, accreting filamentary flows, and represent only a small fraction of the total 
mass inflow rate. The inflowing gas has relatively low metallicity {Z/Zq < 0.2). The outer layers of the 
filamentary infiows are heated due to compression as they approach the disk. In addition to the infiow, 
we find high- velocity, metal-enriched outfiows of hot gas driven by supernova feedback. Our results are 
consistent with observations of halo gas at low z. 

Subject headings: galaxies: evolution - 
methods: numerical 



galaxies: kinematics and dynamics — Galaxy: halo 



1. INTRODUCTION 



face (jHeitsch &: Putman 20091) and the Ha emission along 



Chemical evolution models and analysis of the color 
magnitude diagram of the Hipparcos dataset indicate the 
Milky Way has been forming stars at a nearly constant, 
yet slowly declini ng, rate of 1-3 M(?) yr~^ over the past 



several gigayears (iHernandez et al. 20001 Chiappini et a.1.1 
[2OOII [200a iFuchs et all 120091: iChomiuk fc PovichI I2011D . 
In addition, a continuous supply of low metallicity gas 
coming in at a slightly lower rate is needed to account 
for the metallicity of the long-lived Galactic stars (the G- 
dwarf problem; see lTosilfTQQBt , for a review). The source of 
this fuel has been t he subject of decades of research (see 
iPutman et al.|[2bl2l for a review). 

High-velocity clouds (HVCs) detected in H I surveys 
have long b een suspected as the source of the star for- 
matio n fuel ([Wakker fc van WoerdenlflQQTt iPutman et al.l 
|2003[) . Now that distances to most large complexes are 
known, the mass accretion rate from HVC co mplexes can 
be estimated (IWakkcr 200 1: Wakk er et all [2007, 2008; 
iThom et al.ll2008HPutman et al.ll201^ . The resuhing val- 
ues, however, are too low by about an order of magnitude 
compared to the recent Galactic star formation rate. In 
addition, observations of nearby, star forming spirals of- 
ten show a very limited amount of HI gas in their halos 
(|Sancisi et al.l [200I iHeald et al.l[20TTh . again suggesting 
the HI reservoir in halos may not be the dominant fueling 
source. 

Recently there have been suggestions that extra- 
planar ionized gas may be responsibl e for maintain- 
ing s t ar formation in the Milky Way (jLehner fc Howkl 
1201 It IPutman et al] I2009D . as weh as in other galax- 
ies at intermediate and low redshifts (jBauermeister et al.l 
|2010() . Hydrodynamic simulations that model the mix- 
ing and recooling of cold clouds at the disk-halo inter- 
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the Magellanic Stream ( Bland- Hawthorn et al.|[2007l ) also 
pointed out the potential significance of the influx of warm 
ionized gas. When combining these results with other find- 
ings of substantial amounts of ionized g as in the halos of 
star forming galaxies a t higher redshifts ([Tumlinson et al.l 
l2011l : l^ipp et aLll2011[ ). it is evident that a self-consistent 
dynamical model is necessary to connect galaxy fueling 
mechanisms with the various observational constraints. 

In this paper, we describe a high-resolution cosmologi- 
cal simulation of a Milky Way-mass disk galaxy using an 
adaptive mesh refinement (AMR) code and present its key 
features in terms of the thermal and kinematic distribu- 
tion of gas in such a halo. The high mass resolution {itidm 
and m* ~ lO^M©) and spatial resolution (136-272 pc co- 
moving or better at all times) employed in the simulation 
allow us to study and track the spatial and kinematical 
distribution of the multiphase gas in the halo in great de- 
tail. We describe the simulation in § [5] The results are 
presented in § [31 with the emphasis placed on identifying 
the gas components responsible for inflow onto the galaxy. 
Finally, we examine the evolution of the gas in fllamentary 
flows in the simulation and present a new scenario for gas 
accretion onto Milky Way-sized galaxies in § |4l 

2. METHOD 

We perform simulations with Enzo, an Eulerian hydro- 
dynamics code with AMR capability (Bryan 1999; Nor- 
man fc Bryan 1999; O'Shea et al. 2004). It solves the Eu- 
ler equations using the piecewise-parabolic method (PPM; 
Colella fc Woodward 1984) or the solver used in Zeus 
(Stone fc Norman 1992) to handle compressible flows with 
shocks; we used the latter primarily for numerical stability. 

First, we ran a low-resolution simulation with a peri- 
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odic box of i = 25 /i ^ Mpc comoving on a side with 
cosmological parameters consistent with WMAP5: {flm, 
flA, h, cTg, ris) = (0.279, 0.721, 0.046, 0.70, 0.82, 0.96). 
We identified Local Group-hke volumes by using criteria 
based on the halo mass (mass range 1-2 x 10^^ Mq), the 
mean density (0.60-1.0 times the mean density of the uni- 
verse) and the relatively low velocity dispersion of the ha- 
los (< 200 km s~^) identified within 5 Mpc of a given 
galaxy. We identified four such halos. Then we performed 
a resimulation for one of the four halos using the multi- 
mass initialization technique with four nested levels (five 
including the root grid), achieving to_dm = 1.7 x 10^ Mq, 
within a (~5 Mpc)'^ subvolume. The selected galaxy 
has a halo mass of 1.4 x 10^^ Mq at z = and so contains 
over 8.2 million dark matter particles within the virial ra- 
dius. With a maximum of 10 levels of refinement, the 
maximum spatial resolution stays at 136-272 pc comov- 
ing at all times. Results from the same simulation were 
discussed in Fernandez et al. (2012), particularly in the 
context of H I gas. 

The simulation includes metallicity-dependent cooling 
extended down to 10 K (Dalgarno & McCray 1972), meta- 
galactic UV background, shielding of UV radiation by neu- 
tral hydrogen, and a diffuse form of photoelectric heating 
(Abbott 1982; Joung et al. 2009). The code simultane- 
ously solves a complex chemical network involving multi- 
ple species (e.g., H I, H II, II2, He I, He 11, He III, e^) 
and metal densities explicitly. 

Star formation and stellar feedback, with a minimum 
initial star particle mass of m* = 1.0 x 10^ Mq, are 
also included. Star particles are created in cells that 
satisfy the following two criteria: p > psp and a vio- 
lation of the Truelove criterion (Truelove et al. 1997). 
The star formation efficiency (i.e., the fraction of gaseous 
mass converted to stars per dynamical time) is 0.03 (e.g., 
Krumholz & Tan 2007). Supernovae feedback is mod- 
eled following Cen et al. (2005), with the fraction of 
the stellar rest-mass energy returned to the gas as ther- 
mal energy, esN = 10~^. Feedback energy and ejected 
metals are distributed into 27 local cells centered at the 
star particle in question, weighted by the specific vol- 
ume of the cell. The temporal release of metal-enriched 
gas and thermal energy at time t has the following form: 
f{t,U,U) = il/U)[{t - U)/U]e^p[-{t - where U 

is the formation time of a given star particle, and = 
max(tdyn, 3 x 10^ yr) where idyn = ^Svr/ (32G'/9tot) is the 
dynamical time of the gas from which the star particle 
formed. The metal enrichment inside galaxies and in the 
intergalactic medium (IGM) is followed self-consistently in 
a spatially resolved fashion. 

3. RESULTS 

We extracted a spherical volume from the simulation 
output that extends to the galaxy's virial radius (250 kpc) 
at a uniform spatial resolution of 1.09 kpc/cell. In order to 
examine finer structures, the volume inside 20 kpc in ra- 
dius was extracted at a higher resolution of 0.272 kpc/cell, 
the maximum spatial resolution of the simulation, and this 
replaced the inner volume of the larger sphere. In order to 
focus on gas accretion in the halo, the cylindrical region 
defined by \R\ < 18 kpc and l^;] < 2 kpc whose symme- 
try axis coincides with the rotation axis of the simulated 



disk was removed from this analysis. Hence, the resulting 
quantities reflect the properties of the halo region only. We 
report on our analysis of the simulation result at z = 0, 
unless otherwise specified. The evolution of H I gas in the 
halo at low redshifts ( z < 0.5) was studied in detail by 
[Fernandez et all (|2012l) . 

We find that 70% of the mass influx is concentrated in 
'^17% of the surface area over a large range of radii. This 
implies that the gas inflow occurs along continuous, fila- 
mentary structures. We find three main filaments of warm 
gas that feed the galaxy. Further details on the spatial and 
kinematic properties of these warm filamentary flows will 
be reported in a forthcoming paper. 

3.1. Gas Inflow Velocities 

We examine the distribution of radial velocities of the 
halo gas. The systemic velocity of the galaxy, i.e,, the 
center-of-mass velocity of the dense (n > 0.1 cm^'^) cells 
in the disk, was subtracted from all cells, to focus on the 
relative motion of the halo gas with respect to the galaxy 
itself. 

Figure [1] displays the radial velocity distribution of gas 
in various temperature ranges. It shows the curves repre- 
senting the amount of mass per unit velocity interval in 
three different temperature ranges, plotted against the ra- 
dial velocity. The three temperature ranges were selected 
to be cold (T < 10^ K, blue), representative of H I & Ha 
emission and Lya & Mg II absorbers; warm-hot (10^ K 
< T < 10^ K, yellow), representative of C IV and O VI 
absorbers; and hot (T > 10^ K, red), representative of 
higher level ions such as O VII and O VIII as well as X- 
ray emission. These definitions are used in Figures [2] and 
[3] as well. The gas associated with the last component has 
densities that are usually too low to be detected in current 
observations, except for the region close to the disk. The 
black curve is the sum of the three solid curves mentioned 
above. 

EDITOR: PLACE FIGURE [1] HERE. 

The warm-hot gas dominates the mass over almost the 
entire range of radial velocities. Although the hot gas oc- 
cupies a significant volume fraction, it does not dominate 
the mass because of the low densities. The only exception 
is at the highest radial velocities {vr > 300 km s~^), where 
the hot outfiowing gas contributes ~10* Mq. 

The grey histogram shows the amount of cold gas con- 
tained within 10 kpc of the four gas-rich satellites iden- 
tified within the virial radius of the simulated host halo 
(see below for more details). We picked the radius of 10 
kpc because it is at least 40% (and up to 100%) of the 
virial radii of the satellite subhalos, and so the bulk of the 
cold gas should reside within this volume, unl ess it was 
previ ously ejected or stripped away (see iFernandez et al.l 
120121 ). Three of the four satellites have \vr\ > 100 km s~\ 
suggesting that at least part of the cold gas with extreme 
velocities must be associated with gas contained within or 
stripped recently from the satellite galaxies. 

The mean radial velocity increases with gas tempera- 
ture from more negative to less negative velocities. We 
find that the cold gas has more negative inflow velocities 
(the mass-weighted mean radial velocity (vr) = —82 km 
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s~^) than the warm-hot and hot gases ((fr) = —41 and 
—16 km s~^, respectively). These values are marked by 
vertical lines at the top of the figure. The inflowing veloc- 
ities of cold gas are consistent with observations of HVCs, 
although we leave the details of neutral gas structure, pro- 
jection effects from the position and velocity of the Sun, 
and obscuration by Galactic disk gas to future work. Note 
that the radial velocities alone do not tell us which phase is 
primarily responsible for the gas inflow; we must examine 
the mass flux in the radial direction to answer that. 

3.2. Mass Accretion Rate 

Figure [5] shows the mass accretion rate of gas as a 
function of the galacto-centric distance. To calculate the 
mass accretion rate in thin sphe rical shells centere d on the 
galaxy, we used a formula from lPeek et al.l ()2008[ ): 



^(R.) 



M{R) = 



M, Vi ■ i-h) 
dR 



(1) 



The four sharp features in the H I accretion rate at 
R ~ 63, 78, 188, and 243 kpc correspond to the four gas- 
rich satellites (with Nhi > 10^^ cm~^) found within the 
virial radius at z — 0. The negative values correspond to 
those satellites moving away from the galaxy. The feature 
peaked at i? « 13.5 kpc is also associated with one of the 
four satellite galaxies, although in this case the bulk of its 
mass may come from condensation of gas stripped from the 
satellite galaxy, which is at i? = 78 kpc and moving away 
from the host galaxy at z = ("S19" in [Fernandez et all 
l2012f ). This and other H I features close to the disk may 
indicate cooling of the inflowing halo gas at the disk-halo 
interface. We will investigate this issue in more detail in 
the future. Note that the UV radiation from young stars 
in the disk, which is not included in the simulation, may 
photoionize some of the neutral gas in the halo. Hence, 
our calculated H I mass is an upper limit. 

EDITOR: PLACE FIGURE [3] HERE. 

In Figure [31 each panel corresponds to one of the three 
metallicity ranges: (a) low {Z/Zq < 0.2), (5) intermediate 
(0.2 < Z/Zq < 0.5), and (c) high {Z/Zq > 0.5). The 
black solid curves in the three panels represent the accre- 
tion rate of gas in each of these metallicity bins. To display 
the temperature distribution of the inflowing/outflowing 
gas, the curves are further divided into three temperature 
bins using the same colors as in, e.g., Figure[2Ji- The figure 
demonstrates that it is the low-metallicity gas that dom- 
inates the inflow in all the phases. As the low-metallicity 
gas flows in, it makes a smooth transition from cold to 
warm-hot and then from warm-hot to hot temperatures. 

The bottom panel shows a clear gas outflow of high- 
metallicity gas at the rate of ~1 Mq yr~^. The tem- 
perature of the outflowing gas decreases gradually from 
hot at small radii (R < 100 kpc; bottom panel) to 
warm- hot at larger radii (top right panel), presumably 
due to adiabatic expansion. The metals are carried in 
hot outflowin g gas, as predicted by previous theoreti- 



cal work (e.g., Mac Low et al.lll989l : [Strickland &: StevensI 



[20001: [Marcolini et al.l[200'5r This result is consistent with 



where Mi is the gas mass in the i*^ cell in a given spherical 
shell, Vi is the velocity vector of that cell, fi is the radial 
unit vector, and dR is the thickness of the spherical shell. 
Note that this formula gives the mass accretion rate for 
gas contained in each spherical shell in units of Mq yr~^. 

EDITOR: PLACE FIGURE [2] HERE. 

Plotted in Figure [2] are the net (i.e., inflow minus out- 
flow) mass accretion rates of all gas (a), of the metals 
(6), and of the neutral and ionized hydrogen (c). In each 
panel, the mass accretion rates were divided into the three 
temperature ranges deflned in t j3.1l cold (blue), warm- hot 
(yellow), and hot (red). 

The net mass accretion rate, 3-5 Mq yr^^, at all radii 
is comparable to the star formation rate of the simu- 
lated galaxy at z w (-5 Mq yr-i)I3 The fluctu- 
ation in the mass accretion rate is expected from the 
clumpy and stochastic nature of the accreting mechanisms. 
The amount of neutral gas mass increases at small dis- 
tances from the gala xy (Figure see also Figure 1 in 
[Fernandez et al.[[2O120 . implying some cooling and conden- 
sation of gas close to the disk due to increased background 
pressure. 

The primary result of this paper is shown in Figure 
and^p displaying the mass accretion rates; the overall gas 
accretion is dominated by warm-hot ionized gas, rather 
than cold neutral gas, at almost all radii. The bottom 
panel (Fig. [^t), which shows the accretion rate of hydro- 
gen gas, demonstrates that the ionized gas is responsible 
for most of the mass influx with dMnii/dt « 2-4 Mq 
yr~^, while the neutral gas accounts for only 0.1-0.3 Mq 
yr~^ (excluding satellites). This is due, in part, to the fact 
that the fllamentary flows responsible for roughly half of 
the neutral gas in the halo (jFernandez et al.|[2012[ ) are as- 
sociated with temperatures between 10"* and 10^-^ K and 
are mostly ionized. Closer to the disk (i? < 100 kpc), 
warm-hot gas is gradually heated, and the accretion of hot 
gas becomes increasingly important (but see the caveat in 

As in most other cosmological simulations run to date, the stellar mass (~1.9 X 10" Mq) of our simulated galaxy is too large and too 
centrally concentrated; the associated star formation rate at low redshifts is also too high by factors of 2-3. 



observations finding highly metal-enriched hot gas in the 
X-ray (e.g.. iStrickland fc jie ckman 2007.) and in the ultra- 
violet (e.g.. [Tripp et al.l201ll) . Note that the metals have 
a net outflow rate at almost all radii, although the total 
gas accretion rate always indicates a net inflow. 

4. DISCUSSION 

We showed that the overall gas accretion is dominated 
by warm-hot ionized gas rather than cold neutral gas in 
Milky Way-sized galaxies at low redshifts. According to 
Figure nil, the component that dominates the gas accre- 
tion changes gradually from cold to warm-hot (at R w 240 
kpc) and then from warm-hot to hot (at R ~ 50 kpc), as 
the distance decreases. What is responsible for the gradual 
heating of the inflowing gas? 

To address this question, we plot in Figure [H the cooling 
time and compression time vs. radius for the "hot mode" 
gas (Fig. m upper curves) and "cold mode" gas (lower 
curves), which are hereafter defined as gas with T > 10^-^ 
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K and with T < 1 0^-^ K, resp e ctively . This is similaiB 
to the definition in iKeres et alj (|2005f) . To compute the 
mean timescales, the inverse of the appropriate time was 
weighted by the thermal energy density (|nfcr) in each 
cell and summed over all inflowing (vr < 0) cells in a given 
spherical shell. In computing the cooling time, the diffuse 
photoelectric heating rate was also accounted for. 

For the hot mode inflowing gas, the cooling time (blue) 
is longer than the compression time (red) at all radii, which 
suggests that heating dominates over cooling for this com- 
ponent. On the other hand, if we repeat the same calcu- 
lation for the cold mode gas, the cooling time is shorter 
than the compression time at all radii, leading to net cool- 
ing and condensation of the densest parts of the inflowing 
streams, especially at small radii {R < 100 kpc). If we 
make the same plot for all inflowing gas (bottom panel), 
the two timescales are nearly equal (within a factor of ^2) 
over the entire range in radius, excluding the sites of gas 
associated with the satellites. Overall, the compression 
time is shorter than the cooling time, so heating should 
slightly dominate over cooling. 

EDITOR: PLACE FIGURE H HERE. 

As the gas flows in, its kinetic energy gets slowly con- 
verted to thermal energy due to many weak compressions. 
The energetics work out since Vr ~ 100 km s~^ would 
correspond to ~10^ K in gas temperature. The result is 
consistent with the trend in Figure (TJ i.e., the mean inflow 
velocity decreases as the gas temperature increases. 

Figure [5] shows the distribution of specific entropy, s = 
Tp^~'^ where 7 is the adiabatic index, as a function of 
radius, weighted by gas mass (left panel) and mass flux 
(right panel). (Although the correct term for s is 'adia- 
bat,' we refer to it as entropy following common conven- 
tion.) This shows that the mean entropy is fairly flat with 
radius, increasing by only a factor of two from 20 to 250 
kpc. In the right panel, blue and red represent inflowing 
and outflowing fluxes, respectively. The white contour rep- 
resents zero net mass flux, while the dashed curve shows 
the mass-weighted mean entropy computed from the left 
panel. Comparing the two panels, we see that the inflow- 
ing gas (blue in the right panel) has systematically lower 
entropy than typical gas at that same radius, by a factor 
of 2-5 at i? > 100 kpc. This corresponds to a density 
enhancement of less than a factor of 10 for the inflow- 
ing gas, assuming pressure equilibrium. This is less than 
th e critical cloud ove rdensity required for cooling found 
bv lJoung et al.l (|2012l ) in idealized simulations, and there- 
fore consistent with the fact that we do not sec cooling 
instabilities at large radii {R > 100 kpc). 

At smaller radii, we do see cooling gas - in fact, the 
"beard" in the lower left region of the two plots represent 
the gas cooling at i? < 100 kpc. It is not immediately clear 
how much gas is cooling, but when we rank the individual 
cells by the mass influx at a given distance, we find that 
this cooling gas accounts for only the bottom ~10% of the 



mass influx. The cooling gas corresponds to the innermost 
regions inside the filamentary flows, which have the high- 
est densities and lowest temperatures - hence the lowest 
entropies at a given radius. This gas is responsible for the 
small H I clouds at i? < 100 kpc in Figure 2c. 

The rest of the inflowing gas, with higher specific en- 
tropy, has nearly constant entropy as it flows in; in fact, 
there appears to be a gradual increase in the entropy due 
to mixing or weak shocks. The temperature of this gas 
increases as it approaches the disk and is compressed. It 
might seem surprising that the gas temperature increases 
as the filaments flow in; however, it is consistent with Fig- 
ure m which shows that the cooling time of this inflowing 
warm-hot gas is longer than the compression time for the 
hot mode gas. 

Our resu l t is consistent with previous work (e.g., 
IKeres et all l2005l 2009; Dekcl & Birnboim 200|) that 
found the broad idea that hot mode accretion dominates 
the overall gas accretion in Milky Way-mass halos at low 
z. However, the inflowing filamentary gas is not strongly 
shocked and has lower temperatures than the rest of the 
halo gas. 

EDITOR: PLACE FIGURE [S] HERE. 

The observed neutral gas clouds are only the tip of a 
much larger "iceberg" ; H I clouds represent the small 
densest parts of filamentary fiows that are made up of 
mostly ionized gas, in which they are embedded. Neutral 
gas accounts for a significant but not a dominant frac- 
tion (^1/10) of the accretion rate required to explain the 
current Galactic star formation rate. This is comparable 
to the rece nt estimate of ^0.1 M p, yr^^ from all the HVC 
complexes (jPutman et al.|[20T2l ). Note that the outer enve- 
lope with intermediate velocities will dynamically "shield" 
the H I clouds from the (nearly static) ambient medium 
and decrease the relative velocity, increasing the Kelvin- 
Helmholtz growth time and hence the cloud lifetimes. Fig- 
ure [31 which shows the mass accretion rate broken down 
by temperature for each metallicity range, demonstrates 
that the \ow-Z gas dominates the accretion at all radii for 
all the phases, which supports the filamentary fiow origin 
of the inflowing gas. 

The simulation results are consistent with observations 
of halo gas. The H I clouds found i n galaxy halos are 
largely within 20 kp c of galactic disks (jThom et al.ll20Q8l : 
IWakker et all 120081) . while the (la rgely) ionized c ompo- 
nent extends thro ughout the halo (jProchaska et al.l 1201 it 
iBowen et al.ll2002l ). In addition, the H I component is sur- 
rounded by warm an d warm-hot ga s indic a ting multiphase 
flows are present (jPutman et al.l 120121 : ISembach et al.l 
l2003f ). Finally, substantial quantities of inflowing warm 
gas are consistent w i th th e results of IShull et al.l (|2009D 
and lLehner fc Howkl (l2011h . 

We must point out one caveat in our analysis. The stel- 
lar mass of the simulated galaxy is too concentrated in the 
bulge, and so the gravitational potential well has a slope 

^ In lKeres et al.l l(2005,), the hot mode accretion and cold mode accretion were defined based on the maximum temperature attained by a given 
gas particle. As we cannot follow the history of gas particles in a grid-based code, in order to define the two modes, we use the temperature 
at a given time slice instead of the maximum temperature. For this reason, the mass of the hot mode gas that we compute is a lower limit. In 
particular, some fraction of the gas with T < 10^'^ K at large radii will likely be heated further and so should really count as hot mode gas. On 
the other hand, we believe that the distinction between the two modes at small radii {R < 100 kpc) is robust. Also, note that the definitions 
for the cold mode and hot mode gas are to be distinguished from those for the cold, warm-hot, and hot gas in JjS] 
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that is too steep compared to the Milky Way. For this 
reason, the heating of the incoming gas was hkely overes- 
timated. However, it wih probably introduce only a factor 
of a few error in gas temperature, and we expect the quali- 
tative results reported in this paper to remain unchanged. 

5. CONCLUSIONS 

We analyzed a high-resolution AMR cosmological sim- 
ulation of a Milky Way-mass galaxy including star for- 
mation and supernova feedback, in a fully cosmological 
setting. In summary, our key results are: 

• The inflowing gas is filamentary, and the bulk of 
the infiow is warm-hot (10^ K < T < 10^ K) and 
ionized. 

• Most of the inflowing gas (by mass) does not cool; it 
has nearly constant entropy and so the temperature 
increases as the gas approaches the center. 

• Some of the inflowing gas does manage to cool (in 
the innermost regions of the filaments associated 
with the lowest entropies), but only inside R < 100 
kpc, and mostly within i? < 20 kpc. 



• The inflowing gas has low metallicity [Z/Zq < 0.2). 

• The typical infiow velocities are 50-150 km s~^ and 
generally decrease with increasing gas temperature. 

These results point to a picture in which filamentary gas 
flows, driven by the cosmic web, continue to be important 
in Milky Way-mass galaxies at low redshifts. This inflow 
is not "cold mode accretion" in the sense of iKeres et al.l 
(|200l . since the temperatures typically exceed 10^'^ K 
during the passage through the halo, and radiative cool- 
ing does not dominate heating. However, it also does not 
correspond to classic smooth, hot-mode accretion, and the 
gas in these filaments do not experience a large entropy 
jump at the accretion shock. Instead, this warm-hot fil- 
amentary flow may represent a third mode of accretion 
- important for galaxies like the Milky Way that are not 
far beyond the mass and redshift thresholds below which 
cold-mode accretion dominates. We suggest two areas for 
future work: the fate of the flows as they reach and enter 
the galactic disk at the disk-halo interface, and an explo- 
ration of how mergers, AGN, and feedback may affect the 
gas while it is still in the IGM. 
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Fig. 1. — Mass contained in unit velocity interval as a function of the radial velocity of gas in three different temperature ranges, as 
described in the text, at 2 = 0: cold (T < 10^ K; blue), warm-hot (10^ K < T < 10** K; yellow), and hot (T > 10^ K; red). Positive 
(negative) velocities correspond to outflows (inflows). The black solid line is the sum of the three solid curves in color, while the grey dotted 
line indicates the amount of cold gas contained within 10 kpc of the four gas-rich satellites identified within the virial radius of the simulated 
host halo. Three of the four satellites have \vr\ > 100 km s~^. The warm-hot gas dominates the mass over almost the entire range of radial 
velocities, except at the highest values (i.e., Vr > 300 km s~^), where the hot outflowing gas contributes most significantly with ~10** Mq of 
mass. The mass-weighted mean radial velocities are —82, —41, and —16 km s~^ for the cold, warm-hot, and hot components, respectively, 
and are shown by vertical linos at the top of the plot. 
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Fig. 2. — Net mass accretion rate of (a) all gas, (6) metals, and (c) hydrogen gas as a function of tlie galactocentric distance. The cylindrical 
region containing the disk has been removed from the analysis, as described in the text. In the top and middle panels (a and b), the blue, 
yellow, and red areas indicate the appropriate quantity for the cold, warm-hot, and hot components, respectively. In this figure and Figure [S] 
the sign of the mass flow rate is defined such that infiows will correspond to positive values. In the top panel (a), the dominant component 
for the gas accretion switches from cold to warm-hot gas at ij Ri 240 kpc. Inside of Ri 50 kpc, the hot gas dominates the infiux, although in 
the region closest to the disk {R < 20 kpc), all three components contribute about equally. In contrast to the gas, metals are outfiowing (i.e., 
dMz/dt < in our convention) from the system at almost all radii (6). Despite the low overall mass fraction, the hot component (red), which 
is highly metal-enriched, accounts for the bulk of the metal outfiow. The bottom panel (c) shows that the accretion of the ionized hydrogen 
gas (H II, purple) dominates over that of the neutral hydrogen gas (H I, green). The sharp features at 13.5, 63, 78, 188 and 243 kpc that 
can be identified in all the panels - but mostly distinctively in the H I accretion rate - are from gas associated with the four gas-rich satellite 
galaxies present within the virial radius at ^ = 0. 
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Fig. 3. — Similar to Figure[2]but each panel corresponds to one of the three metallicity ranges: (a) low {Z/Zq < 0.2), (fe) intermediate 
(0.2 < Z/Zq < 0.5), and (c) high {Z/Zq > 0.5). The black solid curves in the three panels represent the accretion rate of gas in each of these 
metallicity bins. To display the temperature distribution of the inflowing/outflowing gas, the curves are further divided into three temperature 
bins using the same colors as in, e.g.. Figure The low-metallicity gas dominates the inflow in all the phases. The bottom panel (c) shows 
a clear gas outflow of high-metallicity gas at the rate of ~1 Mq yr~^. The temperature of the outflowing gas decreases gradually from hot 
at small radii {R < 120 kpc) to warm-hot at larger radii (R > 120 kpc), presumably due to adiabatic expansion. 
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Fig. 4. — (Top panel) Mean cooling time (blue) and compression time (red) as a function of the galactocentric radius, for the inflowing 
hot mode gas (upper curves) and cold mode gas only (lower curves). To compute the mean timescales, the inverse of the appropriate time 
was weighted by the thermal energy density (|nA;T) in each cell and summed over all inflowing (vr < 0) cells in a given spherical shell. For 
the hot mode gas, the mean cooling time is longer than the compression time (tcooi > tcomp, represented by the pink area), leading to net 
heating. For the cold mode gas, on the other hand, the mean cooling time is shorter than the compression time (tcoo! < tcomp, denoted by 
the purple area), resulting in net cooling; note that the difference between the two times increases at small radii. If the same plot is made 
for all inflowing gas (bottom panel), the two timescales are nearly equal (within a factor of ~2) over the entire range in radius. Overall, 
excluding the sites of gas associated with the satellites, the heating should slightly dominate over cooling. The timescales in the top panel 
are reproduced in light grey in the bottom panel. 
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Fig. 5. — Contours showing the distribution of gas mass (left) and the mass fiux (right) in arbitrary units in the specific entropy vs. distance 
plane. The mean entropy is fairly flat with radius, increasing by only a factor of two from 20 to 250 kpc. In the right panel, blue and red 
represent inflowing and outflowing fluxes, respectively. The white contour represents zero net mass flux, while the dashed curve shows the 
mass- weighted mean entropy computed from the left panel. The inflowing gas (blue region in the right panel) has systematically lower entropy 
than typical gas at that same radius, by a factor of 2—5 at i? > 100 kpc, which is less than the critical cloud overdensity required for cooling 
(Joung et al. 2012). Note the "beard" in the lower left corners of both panels; they represent the cooling gas in the innermost regions inside 
the filamentary flows at R < 100 kpc. 



